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We consider the dynamical degenerate four-wave mixing (FWM) model in a cubic nonlinear 
medium including both the time relaxation of the induced nonlinearity and the nonlocal coupling. 
The initial ten-dimensional FWM system can be rewritten as a three-variable intrinsic system 
(namely the intensity pattern, the amplitude of the nonlinearity and the total net gain) which 
is very close to the pumped Maxwell-Bloch system. In the case of a purely nonlocal response the 
initial system reduces to a real damped sine-Gordon (SG) equation. We obtain a new solution of 
this equation in the form of a sech function with a time-dependent coefficient. By applying the 
reductive perturbation method to this damped SG equation, we obtain exactly the cubic complex 
Ginzburg Landau equation (CGL3), but with a time dependence in the loss/gain coefficient. The 
CGL3 describes the properties of the spatially localized interference pattern formed by the FWM. 

PACS numbers: 42.65.-k, 05.45.-a, 89.75.Kd 
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I. INTRODUCTION 

The effect of interaction of light and matter in non- 
linear optics is very often characterized by a coupling 
coefficient which reveals a response of the matter. If 
in addition a mutual mixing of several waves is taken 
into consideration, one deals with nonuniform spatial (or 
spatiotemporal) fields. People usually consider the reac- 
tion of matter is local on the action of the field. But 
this is not always the case. In inertial or nonlocal sys- 
tems the response can be retarded in time or shifted in 
space. As a result, the beam-coupling coefficient takes 
a complex value and some phase addition appears be- 
tween the mixed waves. This can lead to the control of 
parameters of one beam by guiding the properties of an- 
other beam, as well as to the formation of stable localized 
structures (i.e. intensity patterns). In this paper we show 
rigorously that a nonlinear system describing the degen- 
erate wave mixing in a medium which possesses both a 
nonlocal response and relaxation is reduced to one non- 
linear complex Ginzburg-Landau equation (CGLE). We 
develop the technique to obtain the cubic CGLE by us- 
ing the reductive perturbation method for the nonlinear 
dynamical wave coupling system. 

Nonlinear dynamical systems have been studied in- 
tensively during the last decennia after localized struc- 
tures (e.g. solitons) were found in such systems. The 
complex Ginzburg-Landau equation (CGLE) became a 
widely used physical-mathematical model appearing in 
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many branches of physics, chemistry and biology, in or- 
der to describe various localized structures (HH, H, 0]. 
Moreover the CGLE is considered as the simplest model 
containing dissipative soliton solutions, which exist in 
nonequilibrium systems where gain and loss are balanced 
[USUI- I n optics, the dissipative solitons described by 
the Ginzburg-Landau equation appear for pulsed oper- 
ation of passively mode-locked lasers as well as for all- 
optical long-haul soliton transmission lines @, 0, 0, 0] • 

The dissipative models which take into account wave 
interactions have been studied in || first of all as 
the envelope of dissipative solitons emitted by an optical 
parametric oscillator. In [8| the author presents theo- 
retical and experimental studies of stumilated Brillouin 
back-scattering of a continuous pump wave resulting in 
backward-traveling solitary pulses in long fiber-ring cav- 
ities. Nonlinear optical cavities with three-wave inter- 
action in a nonlinear crystal, when the waves have dif- 
ferent frequencies, were considered in Ref. Q. It was 
shown that the spatial dissipative solitons can form spon- 
taneously in that case. For the first time we consider the 
cubic CGLE which appears in the problem of dynami- 
cal interaction of four waves with the same frequencies 
in extended nonlocal media. We show that the CGLE 
is obtained because of a photo-induced nonlocal nonlin- 
ear response which includes a time relaxation term in the 
considered (dissipative) model. 

The next feature that we utilize in the model, the 
nonlocality, reveals itself as an ubiquitous property in 
many branches of physics, e.g. optics, plasmas, Bose- 
Einstein condensates [lCj, [ll( ■ Usually the nonlocal re- 
sponse appears when the nonlinearity is associated with 
some sort of transport process such as heat conduc- 
tion in media with thermal response [l2j], diffusion of 
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molecules or atoms accompanying nonlinear light prop- 
agation in atomic vapors 11 3 , UM . and charge transport 
in photorefractive crystals [l5l Il6j|. Specific properties 
of spatial solitons were investigated in nematic liquid 
crystals, where nonlocal response exists due to reorien- 
tation of anisotropic molecules by a propagating beam 
[l7l [Hi , [T^ . [20l | . The nonlocal nonlinearities with for- 
mation of dissipative optical solitons for a wide-aperture 
laser with saturable absorption were studied recently in 

mm 

One usually investigates stationary changes of the in- 
duced nonlinearity. Our dissipative model includes both 
gain and relaxation of the nonlinearity in a nonlocal 
medium. Since we consider the process of wave coupling, 
the photoinduced nonlocal nonlinearity leads to an effect 
of energy transfer between waves during their propaga- 
tion. In this way the nonlocality plays the role of an 
amplified medium to increase the intensities of one beam 
at the cost of decreasing the energy of another beam. The 
energy transfer effect is observed in the dynamical holog- 
raphy when the interacting beams record a dynamical 
grating, which is shifted from the interference pattern, 
and the same beams diffract from this grating [23l . l24j . 
As a result of this energy transfer both the interference 
pattern and the spatial distribution of the amplitude of 
the nonlinearity get a stable localized pattern along the 
z-longitude direction of the medium [H,[ll,[13- We show 
that the CGLE governs the spatio-temporal dynamics for 
both values. 

Dissipative solitons described by CGLE demonstrate 
a rich variety of unusual properties 0, [28[ , such as sta- 
ble periodic pulsations, bounded solitary waves, periodic 
"explosions" , and collapse. All these unique features may 
find applications in nonlinear wave coupling, in partic- 
ular in the dynamic holography in media with nonlo- 
cal response. Among possible applications in photonics 
let us mention: (i) holographic interferometers including 
phase-shifted interferometers; (ii) traps of light (trapping 
states) in a resonator; (iii) manipulation of pulses hav- 
ing different intensities and durations in order to obtain 
optical logic elements, all-optical switching, pulse retar- 
dation etc., as well as the interaction of pulses not only 
in bulk materials but with thin nonlinear films, nano- 
materials and metamaterials; and many others. During 
the process the medium should possess a nonlocal non- 
linearity, e.g. some kind of transport mechanism; or the 
medium can have a local nonlinearity but a regime of 
moving dynamical gratings should be realized. 

The paper is organized as follows. In section |TT] we 
introduce the four-wave mixing model and recall the ex- 
isting results. In section IIII1 we revisit the derivation of 
the damped sine- Gordon equation and derive a new so- 
lution to the FWM. Finally, in section ITVl we apply the 
method of multiple scale expansion and find as a result 
the cubic CGLE. This procedure proves that the FWM 
as well as the dynamical self-diffraction of waves can be 
considered as a dissipative nonlinear system containing 
stable soliton solutions. 



II. THE INTRISTIC SYSTEM OF THE 
DISSIPATIVE FWM MODEL 

The one-dimensional degenerate FWM initial system 
consists of five partial differential equations, namely, 
four coupled wave equations for slow variable amplitudes 
which connect waves 1 and 2 propagating in a forward 
direction and waves 3 and 4 propagating in a backward 
direction, 

d z A x = -i£A 2 , d z A~2 = i£A 1 , 

d z A 3 = -i£A\, d z A 4 = i£A 3 , (1) 

and the dynamical equation for the medium, which in 
the simplest case includes only a gain being proportional 
to the intensity pattern and an exponential relaxation, in 
the form 

d t £=jl m --, (2) 

T 

We assume here that the interference pattern is formed 
by two pairs of co-propagating waves 

I m = A{A 2 + A 3 A A . (3) 

In Eqs. H])-([3]), Aj(t,z) is the slow variable ampli- 
tude of the j-th plane wave Ej(t, z) = Aj(t, z)e i( - wt ~ S ^, 
£(t, z) is the amplitude of the photoinduced nonlinear 
susceptibility. It must be emphasized that the response 
constant 7 = 7l + *7nl = \lW 9 is complex. The complex 
value of the coupling coefficient £ is an essential feature 
for the existence of soliton-like solutions. The interacting 
waves are connected by the impulse conservation law: 

k\ - k 2 = k A - k 3 . (4) 

We assume the following normalization: all wave ampli- 
tudes are normalized by the square root of the total light 
intensity I = |^i| 2 + \A 2 \ 2 + \A 3 \ 2 + \A 4 \ 2 = const, £ is 
the dimcnsionlcss coefficient of the nonlinearity, z is the 
dimensionless longitudinal coordinate z = (kfj (2k' z ))z' , 
where fco is the amplitude of the wave-vector in the free 
space, z' is the spatial coordinate. We keep the dimen- 
sion of the time-coordinate t in order to display the de- 
pendence of the dispersion relation on the time relaxation 
constant r. In this way, in order to make Eq. ^ dimen- 
sionless, the gain coefficient is normalized by the time 
relaxation constant r and has the dimension [7] = T _1 . 

The system H])-© has been considered for the dy- 
namic holography in the case of a purely nonlocal re- 
sponse 7 = z7nl- Then £ is interpreted as the am- 
plitude of the dynamical grating. As previously found 
[261 I27I [2^ , [30| . the initial system is then reducible to 
a damped sine-Gordon equation (SG), which has a sta- 
tionary solution in the form of a sech function \£\ = 
7C/ cosh [2jCz — p] , with C,p arbitrary constants. Nu- 
merical solutions in the form of periodic oscillations were 
investigated in [30]. The first experimental observation 
of localization of the dynamical grating amplitude along 
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the longitudinal coordinate in bulk ferroelectric crystals 
was made in |27j. For the general case of a complex 7, 
the general stationary solution was later found in [3ll ]. 
together with, in the dynamical case, the general solu- 
tion (expressed with elliptic functions) of the reduction 
(z,t) — > *Jze~ t l T for a purely nonlocal response. 

The ten-dimensional system (HJ-Q is invariant 
under any time-dependent rotation in the space 
{Ai,A 2 , A4, A 3 } which preserves the interference pattern 
©. In a previous work [3l| . we could remove this five- 
parameter unessential freedom and obtain the following 
intrinsic system, 

d z I m = -i£Id, d z I d = -2i£I m + 2i£I m , 

d t £ = jl m , (5) 

T 

admitting the first integral 

4|/ m | 2 +ll = K(t), K arbitrary. (6) 

The real field I d 

I d = -\A 1 \ 2 + \A 2 \ 2 -\A 3 \ 2 + \A 4 \ 2 , (7) 

is the relative net gain. Therefore the four-wave mixing 
is characterized by three intrinsic variables: the intensity 
pattern I m , the grating amplitude £ and the relative net 
gain 7 d . 

This intrinsic system ([5]) is very similar to the pumped 
Maxwcll-Bloch system, an integrable system of nonlinear 
optics defined as [13] 



d x p = Ne, d x p = Ne, 
d x N = -(pe + pe)/2 + 4s, 
dre = p, dre = ~p, 



(8) 



with s a real constant (the system is "pumped" when s 
is nonzero). 

When the four-wave mixing model is undamped (r = 
+00) and has a purely nonlocal response (^(7) = 0), 
while the Maxwell-Bloch system is unpumped (s = 0), 
these two systems can be identified, 



- = 0, 9?(7) = 0, s = : 

T 

z t 2| 7 |/ m 2| 7 |7^ | 7 |/ d -2i£ 2i£ 



X T 



P 



N 



(9) 



and in this case the undamped, purely nonlocal response 
four-wave mixing model admits all the solutions of the 
unpumped complex Maxwell-Bloch system. 



III. DERIVATION OF THE DAMPED 
SINE-GORDON EQUATION 

As shown in [H, [U US H3] , under some specific as- 
sumptions, the system made of the four complex equa- 
tions ([T]) can be integrated explicitly. Because we need it 



later, let us first establish this derivation in full general- 
ity. 

If one represents the complex amplitudes as 
Aj = Mje^i, £ = M e e iip ', {M h M e ,ip h ^) real, (10) 
and introduces the notation 

$ 12 = ip X - Cf 2 - lf e + -, $43 = cp 4 - <f 3 - ip e + — , (11) 



the system ([T]) becomes 



M 2 M e 



Mi 

MiM e 

M 2 
M 3 M e 

Ma 

M 4 M P 



sin $12, 



sin $12, 



sin $ 



(12) 



43 . 



d z M\ = +M 2 M e cos $12, d z ip x 

d z M 2 = -MiM e cos$i2, d z ip 2 

d z M 4 = -M 3 M e cos <Z> 43 , d z <p 4 

d z M 3 = +M 4 M e cos $43, d z ip 3 = + "7^" C sin$ 43 . 
It is then convenient to introduce the first integrals 

fu(t) = l^il 2 + \A 2 \ 2 , fUt) = |^ 4 | 2 + |^ 3 | 2 , (13) 
and to compute the ^-evolution of the two functions 



v 12 = \Atf - \A 2 \ 2 , v A3 = \A 4 \ 2 - |A 3 | 2 . (14) 



One finds 



d z v 12 = +4MiM 2 M e cos$i2, 
d z V4 3 = -4M 4 M 3 M e cos$ 43 , 



(15) 



and, by elimination of Mj, 



(d z v 12 ) 2 = 4 (A 4 2 - v 2 2 ) \£\ 2 cos 2 >S> 12 , 
(d z v 43 ) 2 = 4 (/ 4 4 3 - 4,) \£\ 2 cos 2 $43. 



(16) 



If one defines two functions ui 2 (z , t) , u 43 (z , t) by the re- 
lations 

|£| 2 cos 2 $ 12 - (cU 12 ) 2 , |£| 2 cos 2 $43 = (cU 4 3) 2 , (17) 



the two equations (|16|) can be integrated explicitly in 
terms of the two variables 1*12,^43, 



"12 = -/12 cos(2(wi2 - c 12 )), 

V43 = -/ 2 3 cos ( 2 ( u 43 + C43)), 



(18) 



with C12 and C43 arbitrary functions of t. Basic trigonom- 
etry then yields 



' Ai = +/12 sin(ui2 - c\ 2 )e lVl , 
A 2 = +f 12 cos{u 12 -c 12 y v \ 
M = -/ 4 3 sin(u 43 + c 43 )e llfi , 

k A 3 = +/43 C0S(U 43 + C 4 3)e^ 3 . 



(19) 



We have not succeeded to similarly integrate the equa- 
tions for ifj in (|12p without any additional assumption. 
Let us therefore assume, as was done in (2(| H3, HTJ , that 
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these four equations for the spatial evolution of tpj iden- 
tically vanish, i.e. that sin$ 12 = sin $43 = and the 
phases <pj are independent of z, 

d ziPj = 0, j — 1, 2, 3, 4, 

$12 = n 12 n, $43 = n 43 7r, n 12 and n 43 integers, 
d,tp e = 0, (20) 
and for convenience let us redefine the solution (JTUJ) as 

" £ = (d z u)e i( Pe, 
A 1 = +/i 2 sin(si 2 (M-C 12 ))e^i, 
A 2 = +/i 2 cos( Sl2 ( U -Ci 2 ))e^ 2 , 

A3 = +Us cos(s 43 (u + C 43 ))e^ 3 , 

7T 

$12 = <pi - <fi2 - ¥>e + - = n 12 n 7 S%2 = (-1)" 12 , (21) 

$43 = ¥>4 - ¥3 - + 2 = "43^, «43 = (~1)" 43 , 

Im = \ (A 2 2 sin 2(u - C12) - / 4 2 3 sin 2(u + C43)) 
X e *(^e-7r/2) ; 

^1 = /i 2 2 cos 2(u - C12) + /| 3 cos 2(u + C43), 
, «42,«43 e 2. 

The last complex equation to be enforced ((2]) is equiv- 
alent to the two real equations 

d z d t u + —d z u — K sin(2u + a) = 0, 
r 



7NL smg ^ 2 o _ 2 iC 12 _ f2 o 2iC 43 



12 ,: 



/ 4 V 



(22) 



(9^)(9 t ^ e ) + (cotg 5 )ifsin(2 M + a ) = 0, 7 = | 7 |e 19 . (23) 

If St^e 7^ 0, the ODE (|2"3"|) (with i as a parameter) 
integrates as 

cos(2u + a) = +tanh2 ( K l°° tS 9 (z - z (t)) 



\ ®t<Pe / ...),. 

, „ f K{t) cotgg, , ,.\ 1 J 
sin 2u + a = - sech2 — ^ —(2 - z (i)) , 

then the equation (f22|) restricts this solution to 

tp e = — (t — to): — — KTsech2Kr(z — zq)i(25) 



cos(2m + a) = — tanh2i^T(z — zq), 
sin(2w + a) = — sech2Kr(z — zq), 



(26) 



in which K,to, zo are arbitrary constants. This solution 
can also be viewed as the general solution of the reduction 
I-ax/£ — complex constant of the intrinsic system (0, 



I A = —. 



2K 



= _ e -i(0Otg,)(t-f ]/T jSrT8ech2 ^ _ 

in which the wave number K is arbitrary. Very similar to 
[3lL Eq. (23)], this solution is however new and it depends 
on both space and time. 



■ tanh 2Kt(z — Zq), 



If dt^fe = 0, then 7 must be purely imaginary 

d t ip e = 0, cos <7 = 0, (28) 

this defines the already investigated damped sine-Gordon 
equation. 

The result of the above computation can be summa- 
rized as follows. Under the three assumptions that the 
phases of each Aj are independent of 2, the phase of £ is 
constant, and 7 is purely imaginary, one obtains a solu- 
tion of the system CO)-© represented as (f2Tj) . in terms 
of the real solution u of a damped sine-Gordon equation 
(|2"2"|) (with sin g = 1). The representation (f2"Tj) displays 
the invariance (1, 2, 3, 4, d z , u) — ► (4, 3, 2, 1, — d z , —u) and 
depends on six arbitrary real functions of i (/i 2 , /43, Ci 2 , 
C43J + f2i tfi + V3) and one arbitrary real constant 
(the phase tp e ). 



IV. 5ft( 7 ) = 0. FROM REAL DAMPED 
SINE-GORDON TO CGL3 

It is a classical result 33] that the nonlinear 
Schrodinger equation (NLS) can be derived from the sine- 
Gordon equation by a reductive perturbation method, 
see details in e.g. [3J]. When applied to the real damped 
sine-Gordon equation (f2"2")) , this method yields a complex 
cubic Ginzburg-Landau equation which we now derive. 
Consider the damped sine-Gordon equation (1221) 

E = d t d z u + -d z u - K(t) sin(2u + a(t)) = 0, (29) 

T 

in which u(z,t), K(t), a(t) and r are real. 

Following the classical derivation of NLS from the sine- 
Gordon equation (33l . [33 ] , we define a multiscale expan- 
sion in which u is of order e, while K (t) is of order one, 

, . a(t) 

= e ^2 e 3 <pj{z, ez,..., s k z, . . . , t, . . . , e k t, ...), 

+°° . (30) 
K(t) = s j Kj(sH, . . . , e fe £, . . . ), r = unchanged, 

3=0 
3=0 

and, after renaming the scaled independent variables as 

e k z = Z k , e k t = T k , (31) 

one requires each coefficient Ej to vanish. 
The zero-th order 

Lipo = 0, L = d To d Zo + -8 Zo - 2K (T , ...), (32) 

T 

admits the plane wave-type complex solution 

<p = A(Z 1 ,Z 2 ,T 1 ,T 2 ,...)e i ^ z °- F ( T °>~», (33) 
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in which the complex constant q and the complex func- 
tion F obey the dispersion relation 



i — 

T 



OF 
l dT 



2K = 0. 



(34) 



Since Kq may depend on To, it is convenient to introduce 
the primitive Qq of Kq and to represent the dispersion 
relation by its integrated form 



To Q dQ Q 
t q dT 



(35) 



The physical solution of (|32|) is then chosen as the real 
part of the above complex plane wave 



<p =A(Z x ,Z 2 ,T x ,T 2 ,...)e* +c.c. 

* ry To Qq 

Qq = iq/iQ 2% — . 

t q 



(36) 



The first order equation, which defines the evolution 
of ipi, 



Lip x = -de* - G7e* 
2iK OA 
~q~dZ~ x ' "*8Ti 



Gx = -- 



8A ( 8Qq\ 
l 1^F- - 2 IKi - — ) A, (37) 



requires the vanishing of G x in order to avoid tp x to di- 
verge. This defines two complex conjugate linear PDEs 
for A(Z X ,T X ) and A, and the solution of this first order 
is 



r <pi = o 

A = a(Z x — v g Ti, Z 2 — v g T 2 , T 2 



.Qi — Qo 

-2i , v g 



■2k~n 9Ql '_ (38) 

q 2 ' 9Ti 1: 



in which the complex function of integration a is to be 
determined. 

Since the group velocity v g is generically complex, 
let us introduce the two complex conjugate independent 
variables Xi,Yi, 



X x = Z\ - w g Ti, Y x = X x = Zi - DgTi. 



(39) 



The second order equation similarly defines the evolution 

of (fi 2 , 



Lip 2 



G 2 



^ f e 3(*o+* 1 ) a 3 + e 3(* +* 1 ) s 3 



qG 2 e 



iK x da 2Kq d 2 a 



- qG 2 e 
. da 

'df 2 ~ q 2 dXi ' q 3 d 2 Xx 
4^£ e 2 S? (* +$ 1 )| a |2 a _ £ 



(40) 
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and the cancellation of the secular terms requires G 2 to 
vanish, which defines two complex conjugate nonlinear 
PDEs for a(X u T 2 ) and a(Fi,T 2 ) and yields the value 



<P2 



2tKq 



24tK - 9iq 



-ae 



3(<s> +<s> 



x > + c.c. 



(42) 



Therefore, under the reductive perturbation method, 
the damped sine-Gordon equation (f2"2"| generically yields 
the complex PDE G 2 = Eq. (JH), in which K , K X , K 2 
depend on T 2 . 

In the pure sine-Gordon limit 1/t = 0, K[t) = fco = 
constant, with q real, one checks that the PDE G 2 = 
reduces to the nonlinear Schrodinger equation, 



1 



= 0, K(t) = k , q real 



da 2k d 2 a 



8T 2 q 3 d 2 X x 



■ 4 — \a\ a : 



0. 



(43) 



In the generic case (q complex), the PDE G 2 = 
Eq. (|4"Tj) would be identical to the cubic complex 
Ginzburg-Landau equation (CGL3) if its coefficients were 
independent of T 2 . Let us therefore try to get rid of this 
dependence on T 2 by performing the transformation 



a(X 1 ,T 2 ) = ^,r 1 )e x ^\ 
Z = X 1 -f 1 {T 2 ), r? = / 2 (T 2 ), 



(44) 



in which the complex functions /i , A and the real function 
f 2 can be freely chosen. The best one can achieve is to 
concentrate the dependence on T 2 in only one coefficient, 
e.g. the gain or loss term. Then the functions of the 
transformation are the following, 



A 



2 1 
(I 



K X AT X - 2i 



^ 2 dT 2 . (45) 



The final CGLE is 



T 



M 2 d 2 i> 4 —2—^ ~2<3{q)Z 



1 + 



dr] q 3 d 2 ^ q 



-e r 



+ 2i 



g(g) K 2 
\q\ 2 Kq 



t/> = 0. 



(46) 



in which the coefficient K 2 /Kq depends on T 2 and the 
other coefficients are complex constants. Under the con- 
dition that K 2 /Kq be independent of T 2 , the above PDE 
(|46]) is then identical to the CGL3 equation. 

We want to emphasize here that the CGLE (|46| in- 
cludes only the longitudinal space coordinate £ for the 
variable ip. It does not contain any transverse spatial 
coordinates. 

We thus obtain the CGLE for describing the dynam- 
ics of the FWM in a nonlocal medium with a dissipative 
term, where the dependent variable is the envelope of the 
potential u. With the definition \£\ = d z u and the mul- 
tiscale expansion \£\ = £^2j^o^£j, one can obtain the 
expression which connects the value u with the envelope 
of the spatial distribution of the nonlinearity, 



de 
dZ 

= «#(£,??)e 



(9^0- 



?2 + " 



c.c. (47) 
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Taking into account the equation (|27)) which connects 
£ and the intensity field I m , we obtain that both the non- 
linearity spatial shape and the time behavior of the in- 
tensity pattern are the same spatiotemporal distribution, 
where the magnitudes J m and £ only differ by a constant. 
The complex sign "i" means that the functions I m and £ 
have a relative shift in the spatial coordinate. 

Thus the CGLE (gHD together with the Eqs. (g7J) and 
(j27|) describe the statiotemporal dynamics for the phys- 
ical field I m (the interference pattern), the parameter of 
the phase transition £ and the potential u. 

V. CONCLUSION 

We have obtained the complex Ginzburg-Landau equa- 
tion from the nonlinear systems of the dynamical four- 
wave mixing that includes degenerate wave-coupling in a 
cubic nonlinear medium which has both nonlocal and re- 
laxation response. The obtained CGLE is just the cubic 
one when the response is purely nonlocal, i.e. there is the 
energy transfer only between the interacting waves but 
no phase transfer. In this case the initial FWM system 
is reduced to a damped sine-Gordon equation containing 
the first derivative on the spatial longitudinal coordinate 
z. We show that by applying the reductive perturbation 
method, the real damped sine-Gordon equation reduces 
to the cubic CGLE, except for a loss/gain coefficient de- 
pendent on time. The cubic CGLE describes the dynam- 



ics of the formation of localized states (intensity patterns) 
along longitude z-direction in bulk nonlinear medium. 

The interest is to apply the reductive perturbation 
method to the generic system with the complex response. 
We show the initial generic complex FWM is reduced to 
the intrinsic system, which has three dependent variables 
(I m , £ are complex ones, and Id is real). The intrinsic 
system has a form very similar to the complex Maxwell- 
Bloch system. It coincides completely with the Maxwell- 
Bloch system when at the same time the response is 
purely nonlocal the time relaxation is absent. In op- 
tics there exists an example of reduction of the Maxwell- 
Bloch system to the CGLE, which describes the forma- 
tion of transverse mode structures in lasers [35[ , but they 
are derived in a high-order approximation. 

Till nowadays a number of solutions of the CGLE have 
been found including stable localized patterns Q, i.e. dis- 
sipative solitons. These solutions may find applications 
in the dissipative FWM system. They have great po- 
tential for practical use in photonics by applying wave- 
coupling with a nonlocal medium. 
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